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Abstract 



We study the effect of spatially homogeneous and inhomogeneous shad- 
owing on J/ip production in deuterium-nucleus collisions. We discuss how the 
shadowing and its spatial dependence may be measured by comparing central 
and peripheral dA collisions. These event classes may be selected by using 
gray protons from heavy ion breakup and events where the proton or neutron 
in the deuterium does not interact. We find that inhomogenous shadowing 
has a significant effect on central dA collisions, larger than is observed in 
central AA collisions. The inhomogeneity may be measured by comparing 
the rapidity dependence of J /if) production in central and peripheral colli- 
sions. Results are presented for dAu collisions at ^/snn = 200 GeV and dPb 
collisions at ^/snn = 6.2 TeV. 

I. INTRODUCTION 

The nuclear quark and antiquark distributions have been probed through deep inelastic 
scattering (DIS) of leptons and neutrinos from nuclei. These experiments showed that par- 
ton densities in free protons are modified when bound in the nucleus [1]. This modification, 
referred to collectively as shadowing, depends on the parton momentum fraction x and the 
square of the momentum transfer, Q 2 . Most models of shadowing predict that the modi- 
fication should vary depending on position with in the nucleus [2]. Most DIS experiments 
have been insensitive to this position dependence. However, some spatial inhomogeneity has 
been observed in vN scattering in emulsion [3]. 

Deuterium-nucleus collisions at heavy ion colliders offer a way to measure the structure 
functions of heavy nuclei at higher energies, and hence lower x and higher Q 2 than are 
currently possible with fixed target DIS experiments. In addition, dA collisions are more 
sensitive to the gluon distributions in nuclei than DIS. These dA collisions are preferred 
over pA collisions for technical reasons - the two beams have similar charge Z to mass A 
ratios, simplifying the magnetic optics. Several groups have previously considered structure 
function measurements using pA and dA collisions. However, the spatial dependence has 
not yet been explored. A measurement of the spatial dependence would be of great interest 
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both as a probe of the shadowing mechanism and as an important input to studies of AA 
collisions. 

In this letter, we show that dA collisions can be used to study the spatial dependence 
of nuclear gluon shadowing. We consider two concrete examples: dAu interactions at 
v /5^7 = 200 GeV at RHIC and 6.2 TeV dVb interactions at the LHC. The RHIC con- 
ditions correspond to the Spring, 2003 run. 

We choose the J/i[) as a specific example since it has already been observed in pp and 
Au+Au interactions at RHIC [4]. Since we wish to focus on the magnitude of the shadowing 
effect, we do not include nuclear absorption or transverse momentum, p T , broadening in our 
calculations. The techniques and general conclusions that we discuss should be equally 
applicable to other hard probes such as T, open charm and bottom quarks, Drell-Yan 
dileptons, jets, and dijets. At the LHC, similar approaches should also apply for and 
Z° production. By using a variety of probes, both quark and gluon shadowing can be 
measured over a wide range of x and Q 2 . 

Our J/ip calculations employ the color evaporation model (CEM) which treats all char- 
monium production identically to cc production below the DD threshold, neglecting the 
color and spin of the produced cc pair. The leading order (LO) rapidity distributions of 
J/V's produced in dA collisions at impact parameter b is 



da 
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The LO partonic cc cross sections are defined in Ref. [5], M 2 = xix 2 s and xi^ = 
[Mj \/snn) exp(±|/) ss (mj/^/y/SNN) exp(±y) where mj/^ is the J/ip mass. The fraction of 
cc pairs below the DD threshold that become J/^'s, Fj/^, is fixed at next-to- leading order 
(NLO) [6]. Both this fraction and the theoretical K factor, K th , drop out of the ratios. We 
use m c = 1.2 GeV and Q = M [6]. 

We assume that the nuclear parton densities, F t A , are the product of the nucleon den- 
sity in the nucleus Pa(s), the nucleon parton density f^(x,Q 2 ), and a shadowing function 
Sp S (A, x, Q 2 , r, z) where r and z are the transverse and longitudinal location of the parton 

in position space with s = \/r 2 + z 2 and s' = \J\b — f| 2 + z' 2 . The first subscript, P, refers 
to the choice of shadowing parameterization, while the second, S, refers to the spatial de- 
pendence. All available shadowing parameterizations ignore effects in deuterium. However, 
we take the proton and neutron numbers of both nuclei into account. Thus, 

F t d (x,Q 2 ,f,z)= Pd (s)f t N (x,Q 2 ) (2) 
Ffix, Q 2 , b - f, z') = P a(s')S^ s (A, x, Q 2 ,b- f, z')ff (x, Q 2 ) . (3) 

In the absence of nuclear modifications, Sp S (A, x, Q 2 , r, z) = 1. The nucleon densities of 
the heavy nucleus are assumed to be Woods-Saxon distributions based on measurements of 
the nuclear charge distributions with R Au = 6.38 fm and R Ph = 6.62 fm [7]. The deuteron 
wave function is predominantly S wave with a small D wave admixture. We use the Hulthen 
wave function [8] to calculate the deuteron density distribution. The central densities are 
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normalized so that / d 2 rdzpA(s) = A. We employ the MRST LO parton densities [9] for the 
free nucleon. 

We have chosen two recent parameterizations of the nuclear shadowing effect which 
cover extremes of gluon shadowing at low x. The Eskola et at. parameterization, EKS98, 
is based on the GRV LO [10] parton densities. Valence quark shadowing is identical for 
u and d quarks. Likewise, the shadowing of u and d quarks are identical. Shadowing of 
the heavier flavor sea, s and higher, is calculated separately. The shadowing ratios for 
each parton type are evolved to LO for 1.5 < Q < 100 GeV and are valid for x > 10~ 6 
[11,12]. Interpolation in nuclear mass number allows results to be obtained for any input 
A. The parameterization by Frankfurt, Guzey and Strikman, denoted FGS here, combines 
Gribov theory with hard diffraction [13]. It is based on the CTEQ5M [14] parton densities 
and evolves each parton species separately to NLO for 2 < Q < 100 GeV. Although the 
given x range is 10~ 5 < x < 0.95, the sea quark and gluon ratios are unity for x > 0.2. 
The EKS98 valence quark shadowing ratios are used as input since Gribov theory does not 
predict valence shadowing. The parameterization is available for four different values of A: 
16, 40, 110 and 206. We therefore use the FGS parameterization for A = 206 with both Au 
and Pb. 

Figure 1 compares the homogeneous ratios, Seks98 an d Sfgs f° r Q — 2m c . The FGS 
calculation predicts far more shadowing at small x. The difference is especially large for 
gluons. The FGS approach thus predicts larger antishadowing at x ~ 0.1. The valence and 
sea ratios are also shown, for completeness. 

We have calculated homogeneous shadowing effects on J /if) production at NLO. How- 
ever, the additional integrals required to calculate inhomogeneous effects at NLO lead to 
significant numerical difficulties. Therefore we have compared the LO and NLO calcula- 
tions of homogeneous shadowing in pAu interactions at y/SNN = 200 GeV using the EKS98 
parameterization to verify that the results are essentially identical. 

Some differences may arise from the use of the MRST HO [15] parton densities in the 
NLO code while the LO calculation uses the MRST LO [9] parton densities. In addition, at 
NLO, a new production channel, qg scattering, also contributes. However, its contribution 
is small. The NLO and LO ratios agree well because J/ip production is dominated by gluons 
for xp < 0.6. Thus the gg channel is dominant for \y\ < 3 at RHIC and \y\ < 6.4 at the LHC 
and the LO calculation should reproduce the essential effects of inhomogeneous shadowing. 

Shadowing is not the only nuclear effect on J/ip production. Broadening of the pt 
distributions, observed in fixed-target measurements [16], arises from multiple scattering of 
the initial-state partons [17] and does not affect the total cross section. Thus, a p T -integrated 
measurement should be insensitive to this effect. Nucleon absorption depends on the J/ip 
production mechanism [18]. Absorption is generally expected to increase at large negative 
rapidities where the J/ip may hadronize inside the target. However, this regime is reached 
outside the rapidity coverage of the present detector configurations. 

We now turn to the spatial dependence of the shadowing. Since the data are limited, we 
use parameterizations of the spatial dependence. We compare two parameterizations that 
we have used for inhomogeneous shadowing in AA collisions [2,19-21]. The first, Sp,wSj 
assumes that shadowing is proportional to the local density, p(s), 

S* P , ws (A,x,Q 2 ,r,z) = l + N ws [Sl(A,x,Q 2 ) - 1]^ , (4) 

Po 
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where po is the central density and iVws is chosen so that (I /A) J d 2 rdzpA(s)S PWS = S P . 
When r ^> Ra, the nucleons behave as free particles while in the center of the nucleus, the 
modifications are larger than the average value S P . 

If, on the other hand, shadowing stems from multiple interactions of the incident par- 
ton [22], parton-parton interactions are spread longitudinally over the coherence length, 
l c = Xjlm^x, where is the nucleon mass [23]. For x < 0.016, l c > Ra for any A. The 
interaction is then delocalized over the entire trajectory so that the incident parton inter- 
acts coherently with all the target partons along its path length. At large x, l c <C Ra and 
shadowing is proportional to the local density, as above. In practice, the numerical differ- 
ences between the two models are small and the available data [3] is inadequate to tell the 
difference. However, the both models predict significant differences with the homogenous 
shadowing that is usually used. 

Because of the difficulty of matching the shadowing at large and small x while conserving 
baryon number and momentum, we consider the small x and large x limits separately. This 
matching is necessary because, for y < — 2 at RHIC, l c < Ra u - This regime is measurable in 
the PHENIX muon spectrometers [4]. Equation (4) corresponds to the large x limit. When 
l c ^> Ra, the spatial dependence may be parameterized as 

S P /A, x, Q 2 , f, z) = l + N p (S P (A, x, Q 2 ) - l) {f PA % Z \ ■ (5) 

The integral over z includes the material traversed by the incident nucleon. The normaliza- 
tion requires (1/A) J d 2 rdzp A (s)S P p = S l P . We find N p > N ws . 

Figures 2 and 3 show the ratios of J/ip production per nucleon in dA interactions with 
and without shadowing as a function of rapidity at RHIC {^/snn = 200 GeV) and the LHC 
(\/snn = 6.2 TeV). These ratios are essentially equivalent to the dA per nucleon to pp ratio 
at the same energy since isospin has a negligible effect on J/ip production. 

We consider three cases: central collisions, b < 0.2Ra, peripheral collisions, 0.9Ra < 
b < 1AR A , and minimum-bias interactions, covering all b. Central collisions can be selected 
by choosing events with a large number of 'grey protons' which move with approximately 
the beam rapidity. Grey protons, slow in the target rest frame, are ejected from the heavy 
nucleus in the collision. The name originates from the appearance of their tracks in nuclear 
emulsion. Their momentum must be larger than the Fermi momentum of the target. The 
number of grey protons is related to the number of individual NN collisions, v, which is 
closely related to the impact parameter. Thus the number of grey protons could be used to 
select quite central collisions. This technique was used in the vA study of inhomogeneity 
[3] as well as in lower energy pA studies [24]. In the multiple interaction picture, shadowing 
should be directly proportional to v. 

At a heavy ion collider, grey protons can be detected in a forward proton calorimeter 
(FPC). Neutrons detected in a zero degree calorimeter (ZDC), could also be used, although, 
because of the transverse momentum, the acceptance would be limited. 

The number of iViV collisions is larger in dA than in pA interactions, reducing the 
uncertainty in the impact parameter determination. In a large acceptance detector, it might 
also be possible to use total multiplicity to enhance the impact parameter determination. 

Peripheral collisions could be selected using two criteria. Selecting events with a small 
number of grey protons preferentially chooses a small number of interactions. It is also 
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possible to select events where only one of the nucleons in the deuteron interacts while the 
other is unaffected. Non-interacting neutrons can be detected in a ZDC while non-interacting 
protons could be seen in a FPC. In these NA collisions, the interacting deuterium nucleon 
is usually near the surface of the heavy ion, so that b > R A — 2 fm. However, since the 
deuteron has a significant density even at large distances, not all single-nucleon interactions 
are very peripheral. 

A complication for dA collisions is that the difference between one or two deuterium 
participants in the interaction is significant. This complication can be partly controlled by 
detecting the non-interacting nucleon in forward calorimeters. It can be further controlled by 
studying the J ftp rapidity distribution. The rapidity, y = y c , at which shadowing disappears, 
Sp — > 1, could be used as a calibration point. If the deuterium comes from negative rapidity, 
J/ijj production at negative y occurs via small x partons in the deuteron. The rapidity 
distribution ratios are essentially mirror images of the shadowing parameterizations in Fig. 1. 
Since Sp « 1 at x « 0.025, y c « -0.2 at RHIC and -3.2 at the LHC. 

The J / ip rate at y = y c could be used to normalize the cross sections in different impact 
parameter classes, necessary if the absolute efficiency of cuts on grey protons is difficult to 
assess. Note that y > y c corresponds to low x in the heavy nucleus while at y < y c , the 
antishadowing region is entered. The point y = y c could also be used to normalize the 
dA data to pp data at the same energy, allowing a direct measurement of shadowing since 
higher order corrections unrelated to shadowing cancel out in the ratio. Although y c may 
not be perfectly known, the exact value of y c is not important since a slight error will only 
affect the absolute values of the ratio. Another advantage to this normalization is that, if, 
as expected, the J/ip absorption is independent of rapidity in the region of interest, any 
impact parameter dependence in the absorption correction will drop out. 

The homogeneous results are indicated by the solid (EKS98) and dashed (FGS) lines. 
The local density approximation of inhomogeneous shadowing is shown by the circles and 
squares while the path length approximation is shown by the x's and diamonds. It is obvious 
that Sp has a stronger spatial dependence. 

In central collisions, the inhomogeneous shadowing is stronger. The stronger the homo- 
geneous shadowing, the larger the effect of the spatial dependence. In peripheral collisions, 
the inhomogenous effects are somewhat weaker than the homogenous results but some shad- 
owing is still present. Shadowing persists in part because the density in a heavy nucleus is 
large and approximately constant except close to the surface and partly because the deuteron 
wave function has a non-zero amplitude at quite large distances. Due to these tails, we ex- 
pect some shadowing effects even at very large impact parameters where there are almost 
no interactions. Figures 2(c) and 3(c) show that the impact parameter integrated results 
are identical to the homogenous results, as expected. 

In dA collisions, the rapidity-dependent inhomogeneity is somewhat stronger than that 
predicted for AA collisions [2] . Central dA collisions are confined to the center of the target 
nucleus, where the effect is stronger, while central AA collisions include both central and 
peripheral nucleon-nucleon interactions. 

The predicted spatial dependence would be even stronger for proton projectiles than for 
deuterons since the proton is effectively a point particle in this approach. In contrast, the 
deuteron is the same size, r m 2 fm, as an intermediate mass nucleus. However, it would 
be considerably more difficult to determine the impact parameter in pA collisions since the 



5 



number of grey protons is reduced and no spectator nucleon would be left on the proton side. 
In fact, light ion projectiles may offer an attractive combination, with enough interactions 
for an accurate impact parameter measurement, but little projectile shadowing. 

It should be possible to measure shadowing at midrapidity at RHIC, particularly if the 
gluon shadowing is strong, as with the FGS parameterization. However, measurements at 
high \y\ could prove more insightful, as seen in Fig. 2. At RHIC, measurements have been 
made in pp, dAu and Au+Au at the same y / J/v r /v, allowing direct comparison. The STAR 
electromagnetic calorimeter [25] covers —2 < rj < 1, encompassing part of the antishadowing 
and shadowing regions. The PHENIX muon arms [4] cover 1 < \y\ < 2.4 and could measure 
the low and high x regions while the PHENIX electron spectrometers, \y\ < 0.35, cover the 
region near y c . It would then be possible to compare the antishadowing region, y ~ —2 with 
the shadowing region, y ~ 2. Both detectors should be able to study most of the J/ip pt 
spectrum. A cut on minimum lepton px might reduce the sensitivity to low px J/ips. A cut 
on J/tp pt would increase the x and Q 2 values probed at a given y, allowing cross-checks on 
the shadowing data. The 'cost' of this cut would be some sensitivity to px broadening. 

The higher energy of the LHC moves the antishadowing region out of the reach of any 
planned lepton coverage. Only the low x and transition regions might be observable at 
low pr, see Fig. 3. Both the central detectors of ALICE [26] (|y| < 0.9) and CMS [27] 
(\y\ < 2.4 for the central and forward arms) should be able to clearly distinguish both the 
homogeneous effect and the gross features of the inhomogeneity. The ALICE forward muon 
spectrometer [28], 2.5 < y < 4, could probe the lowest and highest x values possible at the 
LHC, if both dA and Ad collisions were studied. A px cut on the J/ifj would increase the 
x and Q 2 values, perhaps moving the antishadowing region within reach. The T has a very 
similar production mechanism, but would probe larger x and Q 2 values, and should thus be 
sensitive to antishadowing at large \y\. 

In conclusion, we have shown how nuclear shadowing can be measured with J/ip pro- 
duction in dA collisions. The J/ip rapidity is related to the momentum fraction of the 
target parton, allowing for measurement of the x dependence of shadowing. By selecting 
events with different numbers of grey protons and/or a surviving spectator nucleon from the 
deuteron, it is possible to study the spatial dependence of shadowing. In all of the models 
considered, the spatial dependence is large enough to be easily measured. Although we have 
focused on J/i/j, the approach used here should apply to other hard probes such as T, W ± , 
Z°, Drell-Yan, open heavy flavor and jet production. 

We thank K.J. Eskola and V. Guzey for providing the shadowing routines and for dis- 
cussions. This work was supported in part by the Division of Nuclear Physics of the Office 
of High Energy and Nuclear Physics of the U. S. Department of Energy under Contract 
Number DE-AC03-76SF00098. 
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FIG. 1. The EKS98 and FGS shadowing parameterizations are compared at the scale 
H = 2m c = 2.4 GeV. The solid curves are the EKS98 parameterization, the dashed, FGS. 
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FIG. 2. The J/ip shadowing ratio at 200 GeV as a function of rapidity. The results are shown 
for the EKS98 (solid line for homogeneous shadowing, circles and x's for Seks98 ws anc ^ *^EKS98 p 
respectively) and FGS (dashed line for homogeneous shadowing, squares and diamonds for Sp GS ws 
and <Sfgs,p respectively). The impact parameter bins are (a) b/R,A < 0.2, (b) 0.9 < b/2B,A < 1.1, 
and (c) all b. 
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FIG. 3. The J/tp shadowing ratio at 6.2 TeV as a function of rapidity. The results are shown 
for the EKS98 (solid line for homogeneous shadowing, circles and x's for •S'eksqs ws anc ^ *^EKS98 p 
respectively) and FGS (dashed line for homogeneous shadowing, squares and diamonds for Sp G § ws 
and <Sfgs,p respectively). The impact parameter bins are (a) b/RA < 0.2, (b) 0.9 < b/2R,A < 1-1, 
and (c) all b. 
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